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基础 研究 


基于 全 变 分 正则 化 去 爷 积 的 PET 部 分 容积 校正 
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摘要 :目的 提出 全 变 分 正则 化 去 卷 积 方法 并 应 用 于 PET 图 像 前 


分 容积 校正 。 方 法 将 全 变 分 (toal variation, TV) 引 入 到 图 像 退 
化 模型 中 ,提出 基于 全 变 分 正则 化 的 Van Cittert(VC) 和 Richardson-Lucy(RL) 去 卷 积 方法 。 结 果 提出 的 方法 分 别 应 用 于 


NCAT 仿 真 数据 NEMA NU4-2008 IQ 物理 体 模 和 肿瘤 小 鼠 数 据 ( 均 采 用 西门 子 小 动物 Inveon PET 扫 描 得 到 )。 相 比 于 传统 VC 
和 RL 去 卷 积 方法 ,本 文 提 出 方法 在 仿真 实验 中 校正 图 像 有 明显 的 去 品 和 保持 边缘 效果 ,同时 在 小 鼠 实 验 中 肿瘤 区 域 的 活 度 值 
增加 率 为 (10+1.8)% 时 ,图 像 标 准 方 差 (standard deviation, SD) 增 加 率 分 别 从 49.98% 下 降 到 14.26% 和 42.76% 下 降 到 4.70%。 结 
论 本 方法 能 够 在 校正 PET 部 分 容积 效应 的 同时 有 效 抑制 噪声 增加 ,可 望 更 为 准确 地 诊断 肿瘤 。 
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Incorporation of TV regularization in deconvolution for partial volume correction in PET 
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Abstract: Objective We propose a method using total variation (TV) regularization in deconvolution for partial volume 
correction in PET imaging. In the degraded image model, we used TV regularization procedure in Van Cittert (VC) and 


Richardson-Lucy (RL) deconvolution algorithms. These methods were tested in simulated NCAT images and images of NEMA 
NU4-2008 IQ phantom and tumor-bearing mouse scanned by Simens Invoen microPET. The simulated experiment and 


tumor-bearing mouse experiment showed that the algorithms using TV regularization provided superior qualitative and 
quantitative appearance compared with traditional VC and RL algorithms. When the mean intensity of the tumor increased by 
(10+1.8)%, the SD increase percentage was decreased from 49.98% to 14.26% and from 42.76% to 4.70%, suggesting the 


efficiency of the proposed algorithms for reducing PVEs in PET. 
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PET 成 像 中 部 分 容积 效应 (partial volume effect， 
PVE) 会 使 图 像 模 糊 , 图 像 中 组 织 活 度 衰减 ,病灶 边界 识 
别 不 清 ,从 而 影响 临床 诊断 "。 因 此 减轻 部 分 容积 效应 
对 PET 图 像 的 影响 ,对 提高 PET 在 肿瘤 诊断 和 指导 治 
疗 上 的 应 用 至 关 重 要 。 当 前 ,诸多 方法 相继 被 提出 以 实 
现 PET 图 像 部 分 容积 效应 校正 (partial volume 
correction, PVC)。 其 中 一 类 方法 是 引入 MR/CT 解 剖 
先 验 对 PET 图 像 进行 部 分 容积 校正 ,如 在 图 像 重建 或 后 
重建 过 程 ,通过 分 割 和 配 准 解剖 图 像 ,引入 解剖 图 像 的 


像 进行 精确 分 割 ,而 解剖 图 像 分 割 尚 无 绝对 和 鲁 棱 的 方 
法 。 此 外 ,此 类 方法 主要 针对 大 脑 区 域 , 极 大 限制 了 此 
类 方法 的 应 用 。 

基于 图 像 后 重建 的 去 卷 积 校正 技术 提供 了 男 一 类 
解决 方法 。 去 卷 积 方法 只 需 考虑 PET 图 像 本 身 信息 , 方 
法 实现 简单 易 行 ,更 重要 的 是 可 以 针对 全 身 区 域 进行 部 
分 容积 校正 ,如 Teo 等 首先 应 用 Van Cittert(VC) 去 卷 
积 算法 在 PET 肿 瘤 图 像 进行 部 分 容积 校正 。 然 而 ,去 
卷 积 过 程 会 引起 高 水 平 噪声 , 故 作者 提出 应 用 这 种 算法 


边缘 或 者 区 域 信息 汪 。 然 而 这 些 方法 首先 要 对 解剖 图 
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只 是 针对 特定 的 感 兴 趣 区 域 (regions of interest， 
ROIs) ,例如 经 过 很 好 勾勒 出 的 肿瘤 区 域 。Tohka 和 
Reihac” 应 用 VC 和 Richardson-Lucy(RL)"* 去 卷 积 算 
法 进行 脑 PET 图 像 的 部 分 容积 校正 ,结果 表明 这 两 种 去 
卷 积 方法 都 会 导致 图 像 噪声 的 增加 。 为 解决 校正 过 程 
中 噪声 增加 的 问题 ,本 文 将 正则 化 项 引入 到 去 卷 积 模型 
中 。 全 变 分 (total variation, TY)" 已 广泛 应 用 于 图 像 处 
理 领 域 , 能 够 在 保持 图 像 边缘 的 同时 取得 恨 好 去 噪 光 

果 。 针 对 PET 部 分 容积 效应 图 像 活 度 衰减 .边缘 模糊 特 
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点 和 传统 去 卷 积 方法 噪声 增加 问题 ,本 研究 将 全 变 分 方 
法 与 传统 迭代 去 卷 积 方法 二 者 有 机 结合 ,提出 基于 全 变 
分 正则 化 的 VC 和 RL 去 卷 积 算法 ,并 与 中 值 完 验 正则 
化 ~ 和 贝 叶 斯 正则 化 "的 去 卷 积 方法 比较 。 


1 方法 
一 般 而 言 ,去 卷 积 方法 的 通用 框架 是 基于 以 下 模 


开 


一 


O(x)=I1(x) ® h(x) + NY) (1) 
式 中 ,0 是 探测 得 到 的 图 像 ,I 是 真实 理想 的 图 像 ,h 
代表 着 点 扩散 函数 (point spread function, PSF) ,NN 为 加 
性 噪声 ,Z 为 空间 坐标 ,@ 是 卷 积 操作 。 
1.1 经 典 去 卷 积 方法 
1.1.1 Van Cittert (VC) 去 卷 积 算法 Van Cittert(VC) 去 
卷 积 算法 是 从 探测 得 到 的 退化 图 像 O 中 壕 代 恢复 真实 
理想 的 图 像 1, 将 (1) 式 改写 成 最 小 二 乘 形 式 为 : 


11.=5 om-a ®nal (2) 
运用 梯度 下 降 法 求解 得 到 : 
TI(x)=7"(x) + a(O(X) -hx) D TN)) (3) 


式 中 ,o 是 收敛 参数 一 般 设 为 1, 产 和 产 分 别 是 第 
n+1T 次 和 第 n 次 迭代 的 结果 ,0O 和 /与 (HT) 式 定义 相同 。 
1.1.2 Richardson-Lucy (RL) 去 卷 积 算法 Richardson 
和 Lucy 提 出 了 RL 去 卷 积 算法 ,依据 贝 叶 斯 理论 ,由 (1) 
式 可 以 得 到 : 

p(D 


pU|O)=p(O|D: DOF (4) 


其 中 PCllO) 是 后 验 概率 ,P(OID) 是 似 然 概率 ,PCD) 是 
目标 图 像 的 先 验 概率 ,p(0) 是 常数 。 

在 去 卷 积 问题 中 ,假设 PET 图 像 统计 过 程 为 泊 松 
过 程 , 故 似 然 概 率 p(O|) 能 写成 : 

O(x) 

p(OID= TI ® DO] 0 8 DW) 05) 

求解 (5) 式 可 以 得 到 需要 最 小 化 求解 的 代价 函数 : 

= 8 DN- OW:log(h B DO] (6) 


通过 (6) 式 最 小 化 求解 ,可 以 得 到 RL 算法 的 最 终 形 


O(W) 

h(x) ® T(x) 

式 中 ,h(-〗) 是 h(X) 的 共 轿 ,7 和 六 分别 是 第 n+1 次 
和 第 n 次 六 代 的 结果 ,OO 入 与 (1) 式 定义 相同 。 

VC 和 RL 两 种 去 卷 积 方法 都 能 用 于 PET 图 像 的 部 
分 容积 校正 ,但 缺点 是 会 引起 高 水 平 的 噪声 。 因 此 ,为 
了 减少 校正 过 程 中 噪声 的 增加 ,需要 在 去 卷 积 步 又 中 引 
入 正则 化 项 。 


P(X)=7(0) ) ® h(-D)] (7) 


1.2 全 变 分 正则 化 

全 变 分 (total variation, TV) 去 噪 模型 是 一 种 应 用 
很 广 的 滤波 方法 ,优点 是 在 保持 边缘 信息 同时 具有 良好 
的 去 噪 效果 ,三 维 情况 下 的 TV 式 为 : 

2 2 2 
TV(D=|IVildxdydz=| 人 再 + $drdyde (8) 

其 与 经 典 去 卷 积 算法 结合 的 正则 化 形式 分 为 两 种 。 

第 一 种 是 与 VC 去 卷 积 算法 结合 ,在 (2) 式 中 引入 
TV 正则 化 项 , 则 变 为 : 

Jm= 5 ow-0® DN + SY) (9) 


运用 梯度 下 降 法 求解 , 则 (9) 式 的 迭代 求解 式 为 : 
六 (= 产 00+a( OO -NX) DB TO) + 
VI"(x) (10) 
ey 
式 中 ,4Aw 是 正则 化 参数 ,4 是 收敛 参数 一 般 设 为 1， 
div 代 表 的 是 散 度 算 子 ,V 为 梯度 算 子 ,其 余 与 (3) 式 定义 
相同 。 
第 二 种 是 与 RL 去 卷 积 算法 结合 ,在 (6) 式 中 引入 
TV 正则 化 项 , 则 可 以 改写 成 : 
Jrm= > (h ® DO -Olog(h ® DO 
+hr > Vi)| 
求解 (11) 式 可 以 得 到 : 
a O(x) 
ms 
T°(x) (12) 
. ,VI"(x) 
1 -hr div( 一 -一 一 ) 
ie] 
式 中 ,各 变量 定义 与 (10) 式 相同 。 
1.3 中 值 先 验 正 则 化 
为 了 抑制 去 卷 积 迭代 过 程 中 噪声 的 增加 ,在 RL 去 
卷 积 框架 下 ,Kirov 提 出 将 中 值 先 验 (median root prior, 
MRP)"“ 应 用 于 去 卷 积 过 程 中 ,用 来 构建 迟 一 步 MRP 算 
法 : 


Qhry div( 


(11) 


@ h(-x)): 


名 +1,dec 从 

三 ( 门 -MI 
] 十 1 

这 里 的 产 -0 是 在 坐标 点 ? 汐 第 n+1 次 去 卷 积 迭 
代 但 没有 进行 MRP 正则 化 的 结果 , 产 和 疡 …( 六 分 别 为 
在 坐标 点 改 M 第 n 次 和 第 n+1 次 去 卷 积 迭 代 进 行 过 MRP 
正则 化 的 结果 。M = Med(7, 声 RR) 是 以 坐标 点 ?为 中 心 ， 
及 为 半径 的 邻 域 的 像素 中 值 , 其 中 [为 先 验 的 权重 。 
1.4 贝 叶 斯 正则 化 

贝 叶 斯 正则 化 去 卷 积 方法 在 Naqa 文 章 中 指 的 是 在 


1"*!(7)= (13) 
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求解 (1) 式 过 程 中 使 用 最 大 后 验方 法 ,引入 自 吴 图 像 作 
为 先 验 ,在 贝 叶 斯 框架 下 ,得 到 贝 叶 斯 去 卷 积 迭 代 公 
式 为 : 


产 +100= 产 CO:exp[( 


O(x) 
h(x) ® T(x) 
式 中 ,各 变量 定义 与 (7) 式 相同 ,公式 中 的 指数 项 强 
制 正 向 收敛 ,从 而 起 到 正则 化 作用 。 
1.5 数值 解法 


-D)@® h(x] (14 


Ce A (Ee 
对 于 三 维 数据 的 的 divCYZC9 ) 计 算 ,我 们 参考 
[vio)| 
Rudin 等 给 出 如 下 数值 解法 : 
VI 
i 
(vie) 
A A + 
A +m(AxTie A Ti) +m(Al Tin, A Ln) 
A77， (15) 
A’—— + 
Oe 十 m(AYL, A Lix) + m(Ai Ti, A La) 


A: 
,CA TY FA TpsA AAA 区 


其 中 ,i=1..N,, j=1...N,, k=1...N. 


图 1 仿真 NCAT 图 像 
Fig.l Simulated NCAT image. A: simulated PET image; B: degraded PET 
image with partial volume effect and noise. 


2.2 物理 体 模 实验 

PET 物 理 体 模 数据 采用 基于 NEMA NU 4-2008 标 
准 的 图 像 质 量 体 模 四。 图 2 为 采用 有 机 玻璃 加 工 的 体 
模 , 其 圆柱 长 度 为 50 mm, 直径 为 30 mm, 包 含 直径 为 
30 mm\ 高 为 30 mm 的 圆柱 形 空 腔 ,还 有 高 20 mm 的 固 
体 部 分 。 其 中 国体 部 分 有 5 个 直径 依次 为 1.2、3、4、 
5 mm 的 可 填充 空 腔 棒 , 空 腔 棱 与 加 柱 形 空 腔 联通 。 

我 们 在 体 模 中 注入 21 ml 活 度 浓度 为 174.8 kBq/ml 
的 下 -FDG 深 液 ,得 到 的 总 活 度 为 3.67 MBq。 在 南方 医 
院 PET 中心 小 动物 Invoen micro PET 上 进行 扫描 。 
PET 系 统 设置 为 默认 设置 ,能 量 和 事件 符合 时 间 分 别 为 


Azlix = hu( 二 Tix + Tor ni), 
AzT =h (FT tlsm), Mil =h(F La + Tygsn) 


l ign b 
i signa + sign bh 


同时 hh,h,,h. 是 指 像素 的 维 数 , 一 般 都 取 1, 在 边界 
点 ,我 们 定义 以 下 关系 : 

Tox = Lon. 主攻 庆 三 ZN ji 

Zoo =Tix, Tiw,r pe = lin (16) 

Lio=Lin, Tw y= Lin. 

计算 平台 是 Intel Core i5-2400 @ 3.10GHz 四 核 处 
理 器 .4GB 内 存 的 PC 机 。Matlab 版 本 为 7.11.0.584 
(R2010b)。 


min(|a 


3 


2 实验 设计 与 评价 
2.1 NCAT 仿 真实 验 

本 图 1A 是 大 小 为 128 x 128 的 NCAT 仿 真 PET 图 
像 ,参数 均 来 自 文献 [13]。 图 1B 是 仿真 存在 部 分 容积 
效应 和 噪声 的 退化 PET 图 像 , 退 化 图 像 的 部 分 容积 效应 
是 在 图 像 域 中 对 体 模 图 像 进 行 半 高 宽 为 6 mm 的 高 斯 
平滑 ,随后 加 入 均值 为 0, 方 差 为 0.1 的 高 斯 噪声 。 


图 2 NEMA NU4-2008 IQ 体 模 
Fig.2 NEMA NU4-2008 image-quality (IQ) 
phantom. 
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350~650 keV 和 3.432 ns。 我 们 对 体 模 进 行 20 min 扫 
摘 ,数据 存储 为 list-mode 格 式 , 然 后 采用 2D 有 序 子 集 
最 大 似 然 法 (2D ordered subset expectation 
maximization, OSEM2D) 重 建 图 像 (16 个 子 集 ,4 次 迭 
代 )。 最 终 得 到 的 三 维 图 像 数 据 是 256 x 256 x 159, 体 素 
大 小 为 0.39 x 0.39 x 0.79 mm。 
2.3 肿瘤 小 筷 实 验 

临床 FDG 数 据 是 鼠 源 EMT6 乳 腺 癌 小 鼠 数 据 。 我 
们 使 用 南方 医院 Invoen microPET 对 注射 "F-FDG 总 活 
度 为 12.95 MBq 的 肿瘤 鼠 进 行 10 min 扫 描 , 扫 描 过 程 
机 器 的 设置 为 默认 设置 。 重 建 算 法 为 OSEM2D ,参数 
设置 为 16 个 子 集 ,4 次 迭代 。 最 终 得 到 的 三 维 图 像 数 据 
是 256 x 256 x 159, 像 素 大 小 为 0.39 x 0.39 x 0.79 mm。 
2.4 定量 评价 指标 

针对 部 分 容积 效应 校正 后 的 图 像 ,分 别 采用 恢复 系 
数 (recovery coefficient, RC) 和 标准 方差 (standard 
deviation, SD ) 进 行 定量 评价 。 
2.4.1 恢复 系数 恢复 系数 定义 为 : 

RC= Xo 

XROr 

式 中 ,Xkoj 为 探测 图 像 感 兴趣 区 域 的 活 度 值 , XRG 
是 真实 图 像 感 兴趣 区 域 的 活 度 值 。 
2.4.2 标准 方差 标准 方差 表示 图 像 的 噪声 水 平 ,定义 
如 下 : 

SD | -2 x 100% (18) 

式 中 ,N 是 感 兴趣 区 域 的 像素 个 数 ,Xj 为 感 兴趣 区 


域 像素 ) (j=1,.…,N ) 处 的 像素 值 ,X 是 感 兴趣 区 域 的 
均值 。 


(17) 


图 3 NCAT 仿 真 图 像 校 正 结果 
Fig.3 Simulated NCAT image results. A: RL; B: VC; C: RL-TV; D: VC-TV. 


3.2 体 模 实验 结果 

图 4 给 出 了 优化 参数 下 不 同方 法 对 体 模 图 像 校 正 
后 的 结果 ,所 有 图 像 均 在 同一 显示 窗 下 显示 。 图 4 中 各 
个 图 像 的 红色 箭头 均 是 指向 最 小 直径 (1 mm) 的 ROI 
( 记 为 ROI 1)。 由 图 4A 中 可 以 看 出 在 未 校正 图 像 的 
ROI 1 在 图 像 中 亮度 值 很 低 ,经 过 不 同方 法 校正 后 可 以 


同时 ,为 了 表示 噪声 的 增加 率 , 定 义 了 标准 方差 的 
增加 率 : 
op 2 SDerea -SDorsna 
De 

式 中 ,下 标 increase 代 表 增 加 率 ,correct 代 表 校 正 
图 像 计算 值 ,original 代 表 校 正 前 图 像 计算 值 。 
2.4.3 参数 优化 针对 重建 得 到 的 3D 体 模 图 像 数据 , 记 
为 Uncorrected data, 分 别 用 六 种 不 同 的 方法 进行 部 分 
容积 校正 。 方 法 1 是 对 图 像 数 据 进 行 VC 去 卷 积 , 记 为 
VC; 方 法 2 是 对 图 像 数 据 进行 RL 去 卷 积 , 记 为 RL; 方 法 
是 引入 TV 正则 化 项 的 VC 去 卷 积 , 记 为 VC-TV; 方 法 
4 是 引入 TV 正则 化 项 的 RL 去 卷 积 , 记 为 RL-TV; 方 法 5 
在 RL 框架 下 引入 中 值 先 验 , 记 为 MRP; 方 法 6 是 对 图 
数据 进行 贝 叶 斯 去 卷 积 , 记 为 Ba。 在 实验 中 ,所 有 方 
的 迭代 步 数 参 数 都 为 20。PSF 作 为 影响 算法 校正 结 
果 的 重要 参数 ,本文 选取 半 高 宽 (full width at half 
maximum, FWHM) 参数 为 1.5 mm" 。RL-TV 和 
VC-TV 的 TV 正则 化 参数 为 0.005;RL-MRP 的 滤波 核 
大 小 设 为 3 x 3 x3,pB 的 取 值 为 0.3。 


x 100% (19) 
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3 结果 
3.1 仿真 实验 结果 

图 3 是 传统 RL 和 VC 去 卷 积 算法 和 引入 TV 正则 
化 项 去 卷 积 算法 的 结果 。 从 图 3A 和 图 3B 可 以 看 出 RL 
和 VC 都 具有 校正 效果 ,但 噪声 水 平 高 于 校正 前 图 像 。 
而 图 3C.D 是 RL-TV 和 VC-TV 结果 ,可 以 看 出 在 保持 
图 像 边 缘 信 息 的 同时 抑制 了 图 像 噪声 的 增加 ,达到 很 好 
的 校正 效果 ,其 中 RL-TV 校正 效果 好 于 VC-TV。 


看 出 ROI 1 处 亮度 值 都 有 提升 。 同 时 ,5 个 ROI 的 图 像 
是 用 RC 来 定量 分 析 , 空 腔 圆 形 图 像 是 用 来 定量 SD。 

图 $ 描 绘 了 优化 参数 后 的 不 同 去 卷 积 方法 在 不 同 
直径 圆 形 感 兴趣 区 域 的 RC 和 SD 增加 率 变化 曲线 。 从 
图 5A 可 以 看 出 ,直径 越 大 的 圆 形 感 兴趣 区 域 的 RC 值 越 
接近 于 1 ,同时 校正 后 的 RC 值 的 增加 率 范 围 为 9.79%+ 
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图 4 不 同 去 卷 积 方法 校正 的 体 模 图 像 


Fig.4 Various phantom images with different PVC algorithm.s The red arrow points to the smallest ROI (1 mm). 
A: Uncorrected data; B: RL; C: VC; D: RL-TV; E: VC-TV; F: Ba; G: MRP 


2.2%(ROI5) 到 (25.9+1.2)%(ROI 1)。 图 5B 展 示 了 不 
同 校正 方法 的 SD 增 加 率 曲 线 , 即 噪声 增加 水 平 。 首 先 ， 
感 兴趣 区 域 的 直径 大 小 与 校正 后 引入 的 噪声 有 较 大 联 
系 , 基 本 上 是 直径 越 大 ,引入 噪声 越 多 ,但 在 ROI 2 处 的 
RL 和 Ba 方法 引入 噪声 反而 较 在 ROI3 和 ROI 4 处 多 。 
其 次 ,针对 不 同 的 去 卷 积 算法 ,VC 的 SD 增加 率 最 大 , 范 
围 是 从 17.4% 到 20.9%。RL 较 VC 表现 更 好 ,噪声 引入 
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较 少 ,但 仍然 保持 较 高 水 平 。 男 一 方面 ,图 5B 中 不 同 的 
正则 化 方法 都 能 在 一 定 程度 抑制 噪声 的 增加 。 其 中 在 
ROI2-ROI5 条件 下 ,MRP 优 于 Ba,TV 正则 化 算法 优 于 
MRP 和 Ba, 而 且 RL-TV 要 优 于 VC-TV。 特 别 地 在 
ROI1 处 ,Ba 优 于 MRP 和 VC-TV。RL-TV 在 RC 相 等 
情况 下 ,噪声 抑制 效果 最 好 ,SD 增加 率 的 范围 为 10.7% 
~13.4%。 


% of SD increase 


MRP VC-TV RL-TV 
Methods 


VC RL Ba 


图 5 不 同 去 卷 积 方法 在 不 同 直径 感 兴趣 区 域 的 RC 和 SD 增加 率 变 化 曲线 
Fig.5 Plots of RC and SD increase percentage using different methods varying ROIs. A is for RC; B is for SD 


increase percentage. 
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3.3 肿瘤 小 所 实验 结果 

图 6 是 不 同 去 卷 积 方法 对 小 鼠 数据 校正 结果 ,所 有 
图 像 都 在 同一 灰 度 窗 下 显示 。 图 6A 中 ROIA 所 在 图 像 
为 矢 状 面 图 (以 下 相同 ),ROIA 为 肿瘤 区 域 ,ROIB 所 在 
图 像 为 横 切 面 图 (以 下 相同 ),ROIB 用 来 测量 噪声 的 增 
加 。 实 验 中 ,我 们 保证 各 个 去 卷 积 校正 方法 对 肿瘤 区 域 


bb 


4 讨论 

最 初 是 在 1979 年 ,Hoffman 等 提出 用 恢复 系数 
来 描述 并 校正 部 分 容积 效应 。 该 法 主要 是 测量 标准 质 
量 体 模 的 不 同 直径 区 域 的 恢复 系数 ,但 是 临床 图 像 的 病 
灶 区 域 形 状 是 不 规则 的 , 体 模 实验 确定 的 恢复 系数 难以 
用 于 临床 。 随 后 有 学 者 将 迭代 去 卷 积 算法 应 用 于 PET 
肿瘤 图 像 的 部 分 容积 校正 。Kirov 等 人 提出 的 基于 RL 


@l»y@ 
A B 


: a 
v © 


加 


D 

F 
图 6 不 同 去 卷 积 方法 校正 的 小 最 图 像 
Fig.6 Corrected mouse tumor images with different 
algorithms. ROI A: sagittal view (tumor)) ROI B: 


transverse view. A: Uncorrected image; B: RL; C: VC; 
D: RL-TV; E: VC-TV; F: Ba; G: MRP. 


即 ROI A 有 近似 一 致 的 活 度 增强 效果 取 (10+1.8)% 增 
加 率 ) ,从 而 考虑 ROIB 中 噪声 增加 情况 。 表 1 给 出 了 
不 同 去 卷 积 方法 对 小 鼠 数 据 校 正 后 具体 活 度 增 加 率 和 
SD 增加 率 的 数值 比较 , 易 见 ,RL-TV 算 法 在 肿瘤 活 度 增 
加 率 最 大 的 条 件 下 ,噪声 增加 率 最 小 。 


去 卷 积 方法 的 MRP 算法 能 够 有 效 校正 网 像 同 时 抑制 品 
声 , 但 是 在 考虑 拓扑 结构 时 需 手动 调节 的 参数 有 9 个 ， 
本 文 只 是 讨论 了 其 中 最 重要 的 一 个 参数 p。 而 本 文 方 
法 ,总 共 只 有 两 个 参数 ,大 大 减少 了 人 为 因素 的 影响 。 
Boussion 等 "提出 在 去 卷 积 过 程 中 引入 小 波 渡 波 , 但 
3D 情 况 下 只 能 对 每 层 图 像 单独 处 理 后 加 和 ,不 能 实现 
纯粹 意义 上 的 3D 人 处理 ,因此 本 文 不 将 此 方法 作为 比 
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表 1 小 鼠 图 像 校 正 后 活 度 和 SD 增 加 率 比 较 
Tab.1 Comparison of intensity and SD increase percentage for 
mouse images 


Methods 加 ROIA ROIB 
Activity Increase (%) SD Increase (%) 
VC 10.88 49.98 
RL 10.47 42.76 
Ba 10.76 35.66 
MRP 11.06 25.39 
VC-TV 8.20 14.26 
RL-TV 11.32 4.70 


较 。Naqa 等 人 提出 Bayesian 正则 化 去 卷 积 方法 ,这 种 
方法 只 是 对 传统 RL 方法 在 形式 上 做 了 修改 ,引入 了 指 
数 项 的 正 向 收敛 ,在 NU4 2008 体 模 实验 和 肿瘤 小 鼠 实 
验 中 都 表明 此 方法 不 能 有 效 去 除 噪 声 。 

本 研究 在 经 典 去 卷 积 方法 的 基础 上 ,引入 全 变 分 正 
则 化 项 ,在 仿真 体 模 实验 结果 图 3 中 视觉 显示 该 法 可 以 
在 保持 图 像 边缘 的 同时 有 效 抑制 噪声 的 增加 。 定 量 实 
验 中 ,全 变 分 正则 化 项 方法 在 校正 效果 相同 (具有 相同 
的 RC 值 或 者 一 致 的 活 度 增 加 率 ) , 相 比 于 其 他 校正 方 
法 ,具有 最 低 的 噪声 增加 率 , 其 中 RL-TV 优 于 VC-TV， 
表明 了 本 文 提出 方法 的 优越 性 。 在 本 文 的 实验 过 程 中 ， 
2 的 选择 是 根据 量化 准则 优化 选取 ,然而 如 何在 临床 应 
用 中 建立 相应 的 优化 参数 选择 ,是 本 文 进一步 需要 开展 
的 研究 。 基 于 此 ,我 们 将 采用 临床 半 定 量 ,如 标准 吸收 
值 或 定量 的 指标 ,提出 自 适应 的 参数 选择 算法 ,更 好 的 
服务 于 临床 。 

综 上 所 述 ,本 文 提出 了 一 种 基于 全 变 分 正则 化 去 卷 
积 的 PET 图 像 部 分 容积 校正 方法 。 针 对 传统 RL 和 VC 
迭代 去 卷 积 方法 引起 的 噪声 增加 问题 ,本 文 提 出 的 TV 
正则 化 去 卷 积 方法 在 仿真 实验 ,物理 体 模 实 验 .肿瘤 小 
鼠 实验 中 都 展示 了 抑制 噪声 增加 ,保持 图 像 边 缘 的 优越 
性 ,尤其 是 RL-TV 方 法 在 文章 中 提 到 的 方法 中 具有 最 
好 的 去 品 效 果 。 
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